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materials  and  utilizing  a  strain  space  formulation  introduced  in  [1Y7  this 
paper  is  concerned  mainly  with  developments  pertaining  to  strain -hardening 
behavior  consisting  of  three  distinct  types  of  material  response,  namely 
hardening,  softening  and  perfectly  plastic  behavior.  It  is  shown  that  such 
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Abstract.  In  the  context  of  a  purely  mechanical,  rate-type  theory  of  elastic- 
plastic  materials  and  utilizing  a  strain  space  formulation  introduced  in  [1] , 
this  paper  is  concerned  mainly  with  developments  pertaining  to  strain-hardening 
behavior  consisting  of  three  distinct  types  of  material  response,  namely 
hardening,  softening  and  perfectly  plastic  behavior.  It  is  shown  that  such 
strain-hardening  behavior  may  be  characterized  by  a  rate-independent  quotient 
of  quantities  occurring  in  the  loading  criteria  of  strain  space  and  the  cor¬ 
responding  loading  conditions  of  stress  space.  With  the  use  of  special  con¬ 
stitutive  equations,  the  predictive  capability  of  the  results  obtained  are 
illustrated  for  strain-hardening  response  and  saturation  hardening  in  a 
uniaxial  tension  test. 
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1. 


Introduction 


Within  the  scope  of  a  rate-type  mechanical  theory  of  elastic-plastic 
materials,  Naghdi  and  Trapp  [1]  have  recently  discussed  the  advantages  of 
formulating  plasticity  theory  relative  to  yield  (or  loading)  surfaces  in 
strain  space  (rather  than  stress  space).  We  adopt  here  the  loading  criteria 
of  the  strain  space  formulation  as  primary  and  derive  the  associated  loading 
conditions  in  stress  space.  By  comparing  the  local  motion  of  the  loading 
surface  in  stress  space  to  that  of  the  loading  surface  in  strain  space  during 
loading,  we  find  that  three  distinct  types  of  material  response  representing 
hardening,  softening  and  perfectly  plastic  behavior  can  be  defined  in  a  natural 
way.  For  convenience,  these  three  types  of  response  will  be  referred  to 
collectively  as  strain-hardening  behavior.  The  development  leading  to  the 
latter,  as  well  as  illustrative  examples  of  the  results  for  special  constitu¬ 
tive  equations,  are  the  main  objectives  of  the  present  paper.  As  in  [l],  we 
confine  attention  to  the  purely  mechanical  theory  of  elastic -plastic  materials, 
and  base  our  development  on  the  rate-type  stress  space  formulation  of  Green 
and  Naghdi  [2,3]^  and  on  the  alternative  strain  space  formulation  introduced 
by  Naghdi  and  Trapp  [1]. 

By  way  of  motivation,  consider  the  response  of  a  typical  ductile  metal  in 
a  one -dimensional  simple  tension  test  in  which  the  strain  may  be  moderately 
large.  Let  e  and  s  stand,  respectively,  for  the  component  e^  of  the  Lagrangian 
strain  tensor  and  the  component  s^  the  syimnetric  Piola-Kirchhoff  stress  tensor. 
Figure  1  shows  a  plot  of  the  stress  s  versus  the  strain  e  for  the  one-dimensional 
homogeneous  simple  tension  test.  From  the  origin  0  to  the  elastic  limit 


*The  theory  proposed  in  [2,3]  is  a  general  thermodynamical  theory  of  elastic- 
plastic  materials.  The  development  in  [1]  is  carried  out  within  a  purely 
mechanical  framework  which  can  readily  be  interpreted  in  terms  of  the 
isothermal  case  of  the  thermodynamical  theory. 


(identified  by  the  point  1)  the  material  is  elastic,  stress  strictly  increases* 
with  strain,  there  is  no  plastic  straining  and  unloading  takes  place  along  1-0.  On 
the  rising  portion  1-3  (excluding  point  3)  of  the  s-e  curve  both  stress  and  plastic 
strain  strictly  increase  with  strain.  Unloading  from  a  point  such  as  2  takes 
place  along  2-2'  leaving  a  plastic  strain  of  amount  02'.  At  point  3>  s  attains 
its  maximum  valued  On  the  falling  portion  3-4-5  (excluding  point  3)  of  the 
s-e  curve,  stress  strictly  decreases  with  strain,  but  plastic  strain  continues 
to  strictly  increase.  Associated  with  each  point  of  the  segment  1-5  in  Fig.  1, 
there  is  a  unique  yield  point  on  the  s-axis  (i.e.,  in  stress  space)  and  a  unique 
yield  point  on  the  e-axis  (i.e.,  in  strain  space).  For  the  points  1,2, 3, 4, 5 
these  are  denoted  by  A^A^A^A^A^.  and  respectively.  The 

points  A^  and  are  the  initial  yield  points.  As  the  segment  1-5  of  the 
stress-strain  curve  is  traversed,  the  locus  of  the  yield  point  on  the  s-axis 
differs  characteristically  from  that  of  the  yield  point  on  the  e-axis ,  in  that 
the  former  reverses  its  direction  of  motion  while  the  latter  does  not. 

The  usual  loading  criteria  of  the  stress  space  formulation  of  plasticity 
theory,  when  applied  to  the  one -dimensional  case  under  discussion,  require  that 
the  plastic  strain  rate  be  nonzero  whenever  the  yield  point  on  the  s-axis  is 
moving  upwards,  and  be  zero  when  it  is  stationary.  It  is  further  stipulated 
that  the  yield  point  on  the  s-axis  cannot  move  downwards  while  tension  is 


being  applied.  These  criteria  are  consistent  with  the  results  of  the  tensile 
test  for  the  rising  portion  1-3  of  the  stress -strain  curve,  both  for  paths  jl 


Recall  that  a  real-valued  function  f  defined  on  some  interval  J  of  the  real 
line  is  increasing  if  f(xo)Sf(x1)  whenever  and  x2  belong  to  &  and  XgSx^. 

A  function  f  is  strictly  increasing  if  f(x2)>f(xi  )  whenever  x2>x1> 

Similarly,  f  is  decreasing  if  x2  5 x^  implies  f (xj) S  f(x^)  and  strictly 
decreasing  if  x2>x^  implies  f(x2)<f(x^). 

^As  was  observed  by  Naghdi  and  Trapp  [1,  p.  789],  the  maximum  of  the  s-e  curve 
corresponds  to  a  point  which  is  still  in  the  rising  portion  of  the  engineering 
stress  (tt)  versus  engineering  strain  (c)  curve.  The  maximum  of  the  tt-«  curve, 
where  necking  begins,  corresponds  to  a  point  on  the  falling  portion  of  the  s-e 
curve. 


the  type  1-2  and  paths  of  the  type  2-2'.  They  also  demand  the  correct  kind  of 
behavior  for  paths  of  the  type  4-4'  issuing  from  points  on  the  falling  portion 
3-5  of  the  stress -strain  curve.  However,  they  are  clearly  inadequate  for  paths 
of  the  type  3-4  because  the  yield  point  on  the  s-axis  does  move  downwards  for 
any  such  path;  and,  as  was  pointed  out  in  [1],  plastic  strain  is  observed  to  be 
strictly  increasing  in  this  region.  On  the  other  hand,  again  with  reference  to 
the  one-dimensional  case  under  discussion,  the  loading  criteria  of  the  strain 
space  formulation  require  that  the  plastic  strain  rate  be  nonzero  whenever  the 
yield  point  on  the  e-axis  is  moving  outwards  and  that  it  be  zero  whenever  this 
yield  point  is  stationary.  It  is  further  required  that  the  yield  point  on  the 
e-axis  cannot  move  inwards  while  extension  is  occurring.  These  requirements  are 
consistent  with  the  behavior  represented  in  Fig.  1.  Thus,  the  plastic  strain 
is  strictly  increasing  along  the  paths  1-2  and  3-4  and  is  constant  along  the 
paths  2-2'  and  4-4'. 

In  order  to  provide  a  background  for  some  subsequent  developments,  it  is 
desirable  to  make  further  observations  regarding  the  stress -strain  curve  in 
Fig.  1.  In  the  context  of  the  classical  infinitesimal  theory,  we  recall  the 
relations 


e  =  e  +  e_ 
e  p 


ee  =  s/E  ’ 


(1) 


6  P 

where  e  and  e  are  abbreviations  for  the  components  e, .  and  e,,  of  the  elastic 
e  p  11  11 

and  plastic  strains,  respectively,  and  E  (>0)  is  Young's  modulus.  We  note 


that 


de  de  .  .  de  . 

de  _  _ e  (  p  de  _  de_  / e,-l 

ds  ~  ds  ds  '  de  ~  ds  '  ds  ' 

e 

de  1 

Now  with  the  use  of  -r-^-=r>0  and  (2)0,  we  have 

as  £i  d 


de 


(2) 


(3) 


Is-  >  0  if  and  only  if  >  0  ,  §£  <  0  if  and  only  if  -p-  <  0 

as  ae  as  ae 

e  e 


On  the  rising  portion  of  the  s-e  curve  ^jj>  0  (or  equivalently  0) ,  on  the 


A  <3  fl  J  A 

falling  portion  — <0  (j£<0)  an<i  jjj  Poirrt  3  becomes  unbounded.  Then,  at  a 


point  A  on  the  s-e  curve,  with  the  help  of  (2)  and  (3)  it  is  readily  seen  that 


de  f>0 
(<0 


if  and  only  if  A  is  on  the  rising  portion  of  the  curve, 
if  and  only  if  A  is  on  the  falling  portion  of  the  curve. 


w 


while  l  +  dep/deg  becomes  unbounded  at  point  3- 


After  recalling  the  main  features  of  the  purely  mechanical  theory  of 


*  A  A 

elastic-plastic  materials  from  [1,2,3]  in  Section  2,  a  quotient  f/g  of 


quantities  which  are  derived  from  the  loading  functions  f  in  stress  space,  and 


g  in  strain  space,  is  introduced.  It  is  noteworthy  that  while  f  involves  the 


time  rate  of  the  stress  tensor  and  g  the  time  rate  of  the  strain  tensor,  the 


A  .A 


quotient  f/g  is  independent  of  rates.  In  the  latter  part  of  Section  2,  using 
an  equation  obtained  with  the  help  of  a  physically  plausible  work  assumption 


introduced  by  Naghdi  and  Trapp  in  [4],  we  derive  a  geometrically  revealing 


A,  A 


expression  for  the  quotient  f/g  [see  Eq.  (32)].  Next  (Section  3)>  in  terms  of 
A,  A 


the  quotient  f/g,  definitions  are  provided  (see  (43))  for  strain- hardening 
behavior,  i.e.,  for  hardening,  softening  and  perfectly  plastic  behavior,  and 
their  geometrical  implications  are  examined.  It  is  demonstrated  that,  while 
during  loading  the  yield  surface  in  strain  space  is  always  moving  outwards 
locally,  the  corresponding  yield  surface  in  stress  space  may  concurrently  be 
moving  outwards,  inwards  or  may  be  stationary  depending  on  whether  the  material 
is  hardening,  softening,  or  exhibiting  perfectly  plastic  behavior.  Because  of 


While  some  of  the  formulas  in  Section  2  may  appear  to  be  repetitions  of  those 
in  [1],  our  starting  point  and  some  of  our  conclusions  differ  from  [1]  and  for 
clarity  we  have  repeated  these  formulas. 


4. 


our  definitions  (43),  a  variety  of  functions  associated  with  material  behavior  and 

A  A  A 

deriving  from  f/g  or  f  are  found  to  be  positive,  negative  or  zero  according  as  a 
material  exhibits  hardening,  softening  or  perfectly  plastic  behavior.  To  avoid 
undue  repetition,  we  introduce  the  abbreviation  (44)  and  denote  such  conditions 
by  the  letter  H.  Any  function  that  satisfies  conditions  H  can  be  used  to  charac¬ 
terize  strain-hardening  behavior.  By  considering  the  limiting  behavior  of  f/g, 
ve  also  examine  (in  the  context  of  the  developments  of  the  present  paper)  the 
phenomenon  of  saturation  hardening  studied  previously  by  Caulk  and  Naghdi  [5]. 
Definitions  for  saturation  behavior  are  given  at  the  end  of  Section  3. 

The  results  in  Sections  2  and  3  hold  in  the  context  of  the  nonlinear  theory, 
but  in  the  remainder  of  the  paper  attention  is  confined  to  small  deformations 
of  elastic-plastic  materials.  In  order  to  demonstrate  the  predictive  capability 
of  the  strain-hardening  characterization  developed  in  Section  3,  special  sets 
of  constitutive  equations  are  utilized  in  Sections  4  and  5  to  discuss,  respec¬ 
tively,  strain-hardening  response  and  saturation  hardening  under  uniaxial 
loading. 

For  the  particular  constitutive  equations  utilized  in  Section  4,  a  rate- 

independent  characterization  of  strain-hardening  behavior  is  provided  in  terms 

of  a  certain  combination  (20+|0)  of  material  constants.  Moreover,  it  is  shown 

that  both  the  time  rate  of  work-hardening  (k)  and  the  time  rate  of  tension  (s) 

A.A 

may  be  used  to  characterize  strain-hardening  behavior.  While  the  quotient  f/g 
involves  the  coefficient  f  as  well  as  the  derivatives  of  strains  with  respect 
to  stress  (see  Eq.  (64)),  it  is  shown  that  for  a  certain  special  case,  the 
quotient  f/g  may  be  expressed  (see  Eq.  (65a))  in  terms  of  quantities  appearing 
in  (2),  (3)  and  (4)  recorded  earlier  in  this  section.  An  examination  of  details 
of  the  solution  in  Section  4  shews  that  in  uniaxial  tension  and  in  the  sense  of 
our  definitions,  linear  elastic  behavior  is  followed  for  perfectly  plastic 
behavior  by  a  horizontal  6tress -strain  curve,  while  hardening  behavior  is 

5- 


represented  by  a  straight  line  lying  above,  and  softening  by  a  straight  line 
lying  below  the  perfectly  plastic  line. 

Finally  in  Section  5  we  consider  another  set  of  constitutive  equations 
having  in  particular  a  loading  function  employed  by  Caulk  and  Naghdi  [5]  in 
their  discussion  of  hardening  response  in  small  deformation  of  metals.  Again 
it  is  shown  that  a  number  of  different  functions  can  be  used  to  characterize 

strain-hardening  behavior.  Moreover,  it  is  demonstrated  that  the  quotient 

A  A  . 

f/g  may  be  calculated  in  uniaxial  tension  from  a  knowledge  of  the  slope  de/ds, 

found  from  the  stress-strain  curve,  and  the  elastic  constants,  namely  Young's 
modulus  E  and  the  shear  modulus  jj.  ,  and  thus  may  be  easily  identified  experi¬ 
mentally.  Although  our  characterization  of  strain-hardening  is,  in  general, 
different  from  that  discussed  previously  by  Caulk  and  Naghdi  [5],  the  two  sets 
of  results  are  in  agreement  for  the  class  of  materials  for  which  detailed 
comparisons  with  experiments  were  undertaken  in  [5].  In  this  connection  see 
Eqs .  (88),  which  also  include  a  simple  expression  in  terms  of  material  coef¬ 
ficients  for  the  saturation  hardening  constant. 


6. 


of  Quantities  Occurring  in  Loading  Criteria 


Let  the  motion  of  a  body  be  referred  to  a  fixed  system  of  rectangular 
Cartesian  axes  and  let  the  position  of  a  typical  particle  in  the  present  con¬ 
figuration  at  time  t  be  designated  by  x^  =  X^(XA,t),  where  is  a  reference 
position  of  the  particle.  Throughout  the  paper,  lower  case  Latin  indices  are 
associated  with  the  spatial  coordinates  x^  and  assume  the  values  1,2,3* 

Similarly,  upper  case  Latin  indices  are  associated  with  the  material  coordinates 
XA  and  take  the  values  1,2,3*  We  also  adopt  the  usual  convention  of  summation 
over  repeated  indices. 

We  define  a  symmetric  Lagrangian  strain  tensor  by  e^  =  where 

FiK =  3*i/aXA  is  "the  deformation  gradient  relative  to  reference  position  and  6^ 

is  the  Kronecker  symbol.  The  six-dimensional  Euclidean  vector  space  formed 

from  the  components  e  is  called  strain  space.  The  components  of  the  symmetric 

KL 

Piola-Kirchhoff  stress  tensor  are  denoted  by  s^  and  the  six-dimensional  Euclidean 

space  formed  from  these  components  is  called  stress  space. 

We  now  summarize  the  main  ingredients  of  the  purely  mechanical  rate-type 

theory  of  a  finitely  deforming  elastic-plastic  solid  and  base  our  treatment  on 

the  work  of  Green  and  Naghdi  [2,3]  and  Naghdi  and  Trapp  [1],  In  addition  to  the 

strain  tensor  e  ,  we  assume  the  existence  of  a  symmetric^  second  order  tensor- 
KL 

valued  function  e?L  =  e£  (X  ,t)  called  the  plastic  strain  at  X.  and  t,  and  a 

KL  KL  A  A 

scalar-valued  function  K  =  K(XA>t)  called  a  measure  of  work-hardening.  It  is 

assumed  that  the  stress  s^  is  given  by  the  constitutive  equation 

SMN  =  SMN^  ’  U  =  teKL,eKL,IC^  » 

and  that  for  fixed  values  of  eL.  and  tC,  (5)n  possesses  an  inverse  of  the  form 

KL  -L 

eMN  =  eMN^  *  ^  =  tSK7.»eta.,K}  •  (6) 


KL  KL 


^In  [4],  Naghdi  and  Trapp  have  shown  that  the  symmetry  of  e^,  follows  from  a 
physically  plausible  work  assumption  which  will  be  discussed  at  the  end  of 
this  section. 


The  response  functions  s^  and  in  (5)  and  (6)  are  taken  to  be  smooth. 

We  admit  the  existence  of  a  continuously  differentiable  scalar-valued 
yield  (or  loading)  function  g(u)  such  that,  for  fixed  values  of  ejj^  and  K, 
the  equation 


g(u)  =  0  (7) 

represents  a  closed  orientable  hypersurface  £6  of  dimension  five 
enclosing  a  region  t  of  strain  space.  The  work-hardening  parameter  tC 
is  chosen  so  that  g(u)<0  for  all  points  in  the  interior  of  the  region  6.  The 
hypersurface  d£  is  called  the  yield  (or  loading)  surface  in  strain  space. 
Corresponding  to  a  motion  X^,  we  may  associate  with  each  particle  of  the 
body  a  continuous  oriented  curve  Cg  in  strain  space.  This  curve  will  be 
called  a  strain  trajectory.  The  strain  trajectories  are  restricted  to  lie 
initially  in  £  or  on  its  surface  SC,  i.e., 


g(u)  5  0 


(8) 


initially  on  Cg. 

•  -n  _  - 

The  constitutive  equations  for  tc  and  e^  are  LI] 


*  '  CKLeKL  * 


(9) 


and 


4KL  =  ( 


0 

if 

g 

<  0 

> 

(a) 

0 

if 

g 

=  0 

A 

and  g  <  0  , 

(b) 

A 

0 

if 

g 

=  0 

and  g  =  0  , 

(c) 

A 

A 

if 

g 

=  0 

and  g  >  0  , 

(d) 

(10) 


where  =  (^(U)  is  a  symmetric  tensor-valued  function,  a  superposed  dot 
indicates  material  time  differentiation, 


8. 


and  where  X=X(U)  and  Pkl=Pkl(^)  are^,  respectively,  a  scalar-valued  function 

A 

and  a  symmetric  tensor-valued  function.  The  quantity  g  is  the  inner  product 

of  the  tangent  vector  e,_.  to  a  strain  trajectory  C  and  the  vector  .  When 

aM  ? 

g  =  0  and  dg/de^^O,  g  gives  the  inner  product  of  e^  and  the  outward 
normal  vector  to  the  yield  surface  96,  where  the  notation  Sg/Se^  stands  for 
the  symmetric  form  1^/96^  +  9g/9eNM) .  The  conditions  involving  g  and  g  in 
(10)  are  the  loading  criteria  of  the  strain  space  formulation.  Using  conven¬ 
tional  terminology,  these  four  conditions  in  the  order  listed  correspond  to 
(a)  an  elastic  state  (or  point  in  strain  space);  (b)  unloading  from  an  elastic- 
plastic  state,  i.e.,  a  point  in  strain  space  for  which  g  =  0;  (c)  neutral 

loading  from  an  elastic-plastic  state;  and  (d)  loading  from  an  elastic-plastic 

A  .  . 

state.  We  assume  that  the  coefficient  of  g  in  (lOd)  is  nonzero  on  the  yield 
surface  and,  without  loss  in  generality,  we  then  set 


Prl  7  0  ,  x  >  o  . 


(12) 


In  order  to  provide  a  geometrical  interpretation  of  the  conditions  (10), 
we  need  to  record  the  material  time  derivative  of  the  loading  function,  namely 

•  <13) 

aeKL 

where  (11)  has  been  used.  It  follows  from  (7),  (9)  and  (10® )  that,  in  an 
elastic  state  the  strain  trajectory  Cg  lies  in  the  interior  of  6  which  is 
referred  to  as  the  elastic  region  in  strain  space,  and  the  yield  surface  96 
remains  stationary.  Similarly,  by  (7),  (9),  (10b)  and  (13) ,  during  uploading 
the  strain  trajectory  intersects  the  yield  surface  96  and  is  moving  in 
an  inwardly  direction,  with  the  function  g  decreasing,  while  96  itself 


^Our  notation  X  corresponds  to  X  in  [1]. 


remains  stationary.  Likewise,  from  (7),  (9),  (10c)  and  (13)  during  neutral 


loading  the  strain  trajectory  Cg  lies  in  the  yield  surface  while  the  latter 
remains  stationary  and  g  =  0.  Finally,  from  (7),  (9),  (lOd)  and  (13)  during 
loading  the  strain  trajectory  Cg  intersects  and  is  moving  in  an  outwardly 
direction.  It  is  stipulated  in  this  case  that  is  locally  pushed  outwards 
by  the  strain  trajectory  Cg  so  that'*' 

k  =  0  ,  (14) 

A 

if  g=0,  g>0.  Thus,  positive  values  of  the  function  g  can  never  be  reached 
on  a  strain  trajectory  and  the  condition  (8)  holds  for  all  time.  It  follows 
from  (9)j  (lOd),  (13)  and  (14)  that  during  loading 

+  ’  0  •  (15) 
aeKL 

Therefore,  since  the  coefficient  of  g  is  independent  of  e^,  we  have 

=  °  (16) 

deKL 

at  all  points  on  the  yield  surface  d£  through  which  loading  can  occur.  We  note 
that  equations  (5)^  and  (6)^  hold  during  loading,  neutral  loading,  unloading 
and  in  an  elastic  state. 

For  a  given  motion  and  an  associated  strain  trajectory  Cg  we  may  utilize 
the  constitutive  equations  (5)^,  (9)  and  (10),  together  with  appropriate  initial 
conditions  for  e^T  and  K,  to  obtain  the  corresponding  stress  trajectory  C  ,  a 
continuous  oriented  curve  in  stress  space.  In  a  similar  fashion  (6)^  may  be  used 
to  obtain  C  from  C  .  Furthermore,  for  a  given  loading  function  g(li),  with  the 

6  5 

aid  of  (6)^,  we  can  obtain  a  corresponding  function  f(lf)  through  the  formula 


In  the  literature  on  plasticity  this  is  called  the  "consistency"  condition, 
namely  that  loading  from  an  elastic-plastic  state  leads  to  another  elastic- 
plastic  state.  For  references  and  background  information  in  the  context  of  a 
stress  space  formulation,  see  for  example  Naghdi  [6,  pages  l4l,  137] • 


10. 

'  i 

‘u  ..  . 


where  the  variables  U  and  V  are  defined  by  (5)g  an<*  (6)g>  respectively. 

Conversely,  (5)-j_may  be  used  to  obtain  g  ftom  f.  Because  of  the  assumed 
smoothness  of  (6). ,  for  fixed  values  of  e£T  and  K,  the  equation 

X  fvL 

f(u)  =  o  (18) 

represents  a  hypersurface  3§  in  stress  space  having  the  same  geometrical 
properties  as  the  hypersurface  36  in  strain  space.  The  region  enclosed  by  3§ 
is  denoted  §.  It  follows  from  (17)  that  a  point  in  strain  space  belongs  to  the 
elastic  region  £  (i.e.,  g(U)<0)  if  and  only  if  the  corresponding  point  in  stress 
space  satisfies  f(v)<0  and  hence  belongs  to  §.  Similarly,  by  (1?)  and  (18)  a 
point  in  strain  space  belongs  to  the  yield  surface  3£  (i.e.,  g  =  0)  if  and  only 

if  the  corresponding  point  in  stress  space  belongs  to  as  (i.e.,  f=0).  Hence,  we 

refer  to  the  interior  of  S  as  the  elastic  region  in  stress  space,  and  to  3?  as  the 

yield  (or  loading)  surface  in  stress  space.  We  have  seen  that  (8)  holds  for 

all  t.  Therefore,  by  (17),  every  stress  trajectory  C  is  restricted 
to  lie  in  S  or  on  its  surface  3S  and  positive  values  of  f  can  never  be  reached. 

We  note  that  any  function  of  variables  U  can  be  written  as  a  different 
function  of  variables  V  and  vice  versa,  e.g.,  =  C^Cu)  =  C^V)  which 

occurs  in  (9). 

In  [1]  a  comparative  basis  was  provided  between  the  two  independent  sets  of 
loading  criteria  for  the  stress  space  and  the  strain  space  formulations.  A 
correspondence  between  the  two  sets  of  loading  criteria  was  established  for  all 
conditions  except  that  during  loading.  The  approach  in  the  present  paper 
differs  from  that  of  [l]  in  that  the  loading  criteria  of  the  strain  space 
formulation  are  regarded  as  primary  and  associated  loading  conditions  in 


stress  space  are  deduced  from  the  former*.  Although  in  the  examination  of 
the  loading  criteria  our  starting  point  and  conclusions  are  different,  the 
arguments  employed  parallel  those  of  [1],  Thus,  taking  the  material  time 
derivative  of  (17)  and  making  use  of  (13)  we  obtain 


f  =  f +  -&L-  4P  l  - 
p  KL  ate 
deKL 


A 

g  +■ 


3e 


:p  +  a& ; 

KL  a<  * 


KL 


where 


(19) 


i  =  -^s 
asMN  m 


(20) 


A  • 

The  quantity  f  is,  of  course,  the  inner  product  of  the  tangent  vector  s^  to  a 

stress  trajectory  Cg  and  the  vector  af/ds^.  In  view  of  (19) 


A  A  •  n  •  .  . 

f  =  g  if  e?L  =  0  and  K  =  0  .  (21) 

rJj 

* 

Considering  an  elastic  state,  f  =  g<0,  (10a)  holds,  K  =  0  by  (9)  and  hence 
A  A  .  . 

f = g  by  (21).  Since  the  yield  surface  3£  in  strain  space  is  stationary  so  also 

is  the  yield  surface  38  in  stress  space.  The  stress  trajectories  remain  in  the 

interior  of  5.  It  is  clear  from  (9),  (10b),  (17)  and  (21)  that*  f  =  0  and  f<0  if 
A  ,  A  . 

g  =  0  and  g<0.  In  this  case  (g  =  0,  g<0)  the  stress  trajectory  C  intersects 

s 

3$  and  is  directed  inwards,  with  the  function  f  decreasing  in  value,  while  38 


*These  derived  conditions  are  not  the  same  as  the  loading  criteria  usually 
assumed  in  the  stress  space  formulation. 

* 

In  [1]  it  was  possible  to  prove  the  converse  of  this  statement  because  of  the 
Independent  loading  criteria  that  were  assumed  in  the  stress-space  formulation. 
It  will  become  clear  presently  that  in  the  context  of  this  paper  the  converse 
statement  does  not  hold. 
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itself  remains  stationary.  It  follows  from  (9),  (10c),  (17)  and  (21)  that  f=0 
and  £=0  if  g  =  0  and  g  =  0.  In  this  case  (g  =  0,  g  =  0)  the  stress  trajectory 
C  lies  in  the  surface  38  which  remains  stationary  and  f=0. 

In  the  case  of  loading  from  an  elastic-plastic  state,  it  follows  from 
(9),  (lOd)  and  (19)  that 

A=1  +  XpKL^  3p  “  p  ^  +  3ic)CKL-*  *  (S  =  0,  g>0)  .  (22) 

«  aeKL  ^KL 

A.A 

In  the  developments  that  follow,  the  quotient  f/g  can  be  expressed  in  a 

number  of  different  forms .  In  order  to  establish  one  such  form  we  note  that 

by  (17),  (5)^}  (6)^  and  the  chain  rule  of  differentiation 

A  A 

ag  af  =  _af_  ^mn  ag  aeMN 


aeKL  aeKL 


deMN 


ag  af  af  °^mn  ag  0  mn 


ate  3k  as^  aic 


^mn  dic 


With  the  use  of  (23),  (22)  can  be  rewritten  as 
£=1+x_  _d£_  r^MN  ,  ^MN  _  , 

J  Xpkl  as  ^  p  aic 
g  MN  ae^ 

•  <e=0’'>0)  •  ( 
KL 

Another  useful  form  of  the  quotient  fi/g  that  may  be  derived  from  (22)  with  the 


(24a) 


help  of  (12)2  and  (16)  is 


C  1  = 

\  n  Wt  * 


A  *pKLl.  p  3K  UKLJ 


0  +  1 
PKL'-.  p  atC  CKLJ 

^KL _ 


>  (g=0,  g>0)  . 


(24b) 


Since  the  right-hand  side  of  (24b)g  is  independent  of  rates,  it  is  clear  that 

the  quotient  f/g  is  independent  of  rates  and  has  the  same  value  for  all  strain 

trajectories  through  a  given  elastic-plastic  point  on  8fi.  Also,  in  view  of  (17), 

A.A 

f/g  is  dimensionless.  Clearly  a  knowledge  of  all  constitutive  equations  is  required 
See  the  previous  footnote. 


.*u- 


A,A 

for  the  calculation  of  f/g. 

We  now  turn  to  the  work  assumption  of  Naghdi  and  Trapp  [4,7]*  Starting 
with  the  assumption  that  the  external  work  done  on  an  elastic-plastic  body  in 
any  smooth  spatially  homogeneous  closed  strain  trajectory  is  nonnegative,  it 
was  demonstrated*  in  [4]  that 


A  A 
dSMN  .  ^MN 


f — M  +  — £2EL  d  IIP  4  SO 
\  p  bK  CKL-’  KL  MN  ~ 
d  KL 

A 

during  loading  or  neutral  loading,  i.e.,  when  g  =  0,  gSO.  In  the  case  of 
neutral  loading  it  follows  from  (10c)  that  the  left-hand  side  of  (25)  vanishes 
and  (25)  is  satisfied  trivially,  while  in  the  case  of  loading,  it  follows  from 


(25) 


(lOd)  and  (12 )g that  (25)  becomes 


A 


A  A 

faSMN  aSMN  _  , 

^  p  dK  CKL^P,"e 
aeKL 


KL  MN 


s  0 


(26) 


with  g  =  0  and  g>0.  The  coefficient  of  e^  in  (26)  is  itself  independent  of 

.  .  A 

e^  and  the  inequality  must  hold  for  all  choices  of  that  satisfy  g>0. 

Therefore,  by  the  same  argument  used  in  Section  5  °f  [4],  we  deduce  that 


(27) 


be 


KL 


'MN 


evaluated  on  the  yield  surface  g  =  0,  where  the  scalar  function  y  satisfies 


Y  =  y  (U)  S  0  . 

We  emphasize  that  (27)  holds  even  for  a  motion  that  is  not  homogeneous*. 


(28) 


In  order  to  compare  (27)  with  the  results  of  [4]  we  multiply  on  both  sides 


of  (27)  by  \g  and  utilize  (10d)  to  obtain 

I-  p  b<  CKLJ  KL  Ky  g 
aeKL 


be 


(g  =  0>  g>0) 


(29) 


MN 


*See  equations  (5*2),  (5-3)  and  (4.11)  of  [4];  the  notation  in  [4]  cor¬ 
responds  to  in  the  present  paper. 

*For  a  discussion  of  this  point,  see  [4,  p.  40]  or  [7,  p.  63]. 


14 


41 

Recalling  the  restriction*  (12^  and  (28),  we  define  a  function  y  by 


*A 

V  =  KY  g  S  0 


and  from  (29)  obtain 

ae^j  MN 

Equation  (31)  is  the  same  as  (5.4)  of  [4].  We  note  that  (31)  involves  rates, 
while  (27)  does  not.  We  have  shown  that  (27)  implies  (31)-  Conversely,  it 
follows  at  once  from  (lOd)  and  (12^ that  (31)  implies  (27)  with  Y*=v/\g,  as 
in  (30). 

From  (27)  and  (24a )^  follows  an  expression  for  f/g  in  the  form: 


|=1-XY*A  >  (g  =  0,  g>0)  ,  (32) 

g 

where 

»  (g  =  f  =  0)  .  (33) 

dMN  deMN 

The  quantity  A  (when  nonzero)  represents  the  inner  product  of  the  normal  to  the  yield 
surface  dS  in  stress  space  and  the  normal  to  the  yield  surface  in  strain  space. 

For  some  purposes  it  is  convenient  to  express  the  constitutive  equation 
(5)^ in  terms  of  an  equivalent  set  of  kinematical  variables  in  the  form 


SMN  =  SMN^eKL'eKL,eKL,K^  '  ^ 

A  ♦ 

Suppose  that  the  partial  derivatives  possess  the  symmetries 

The  function  y  on  the  right-hand  side  of  (31)  depends  on  the  variables  e^, 
ef©I»K  and  eMN"  A 

'This  is  equivalent  to  the  condition  that  be  derivable  ftrcm  a  potential,  as 
indeed  is  the  case  in  the  general  thermodynamical  theory  (see  Section  4  of  (3]) 
of  which  the  present  development  may  be  regarded  as  corresponding  to  the  iso¬ 
thermal  case.  The  existence  of  a  potential  in  the  purely  mechanical  theory 
can  also  be  demonstrated  by  an  argument  based  on  the  work  postulate  of  [4]. 


A  A 
asMW_a,KL 

^KL  J*MN 
obtain* 


,  then  in  a  manner  similar  to  that  in  [4,  Sec.  5]  team  (27)  we 


_2L. 


MN  3e 


KL 


Ri 


,  (g  =  f  =  0)  . 


(35) 


It  is  clear  frem  (35) »  (30),  (lOd)  and  (12)  that  if  the  response  function 
s^  in  (34)  is  independent  of  its  second  and  third  arguments,  i.e.,  if 
3s, 


MN 

5mn 


MN 


(36) 


=°’  lf=0’  then 

36  KL 

PKL-^^0 , .  (6  =  0) 

and  is  directed  along  the  normal  to  the  yield  surface  <bg  in  stress  space, 
as  also  is  during  loading.  It  follows  from  (36)^  and  (28)  that  in  this  case 


Y  >0  ,  (g  =  0) 


(37) 


and  hence,  during  loading,  in  view  of  (30)  and  (12), 


Y  >  0 


(38) 


also.  When  satisfies  (36)^,  (l6)  can  be  written  a* 

1*xv*^r<^"ttfcKL)  -°  -  <«-°>  •  <39) 

^  ^kl 

* 

The  last  result  can  be  used  to  solve  for  the  product  \y  and  (30)  then  gives  y. 

* 

Also,  we  may  set  y  equal  to  an  arbitrary  positive  scalar-valued  function  of  the 
variables  U  and  then  use  (39)  to  determine  \.  Thus,  in  the  special  case  in  which 
in  (34)  depends  only  on  its  first  argument,  no  constitutive  equation  is 
needed  for  p^. 


We 


A,A 

express  f/g  as 


observe  that  when  p„.  satisfies  (36). ,  then  (24b)..  may  be  used  to 

lUi  X  1 


A. A  *  ,  A  . 

f/g  =  4y  r  ,  (g  =  0,  g>0)  , 


*The  symmetry  of  p^  and  hence  ej^  follcws  from  (35).  See  [4,  Sec.  5l» 


(40) 


3.  Strain-Hardening  Response.  Geometrical  Interpretation. 


The  quotient  f/g  which  occurs  in  (24b)2  and  related  equations  in  Section  2, 
is  utilized  here  to  define  three  distinct  types  of  strain-hardening  response  for 
an  elastic-plastic  material.  These  definitions  are  as  follows:  An  elastic-plastic 

material  is  said  to  be  hardening,  softening  or  exhibiting  perfectly  plastic 

,  A  * 

behavior  during  loading  (g  =  0,  g>0)  according  to  whether 

(a)  f/g  >  0  (for  hardening)  , 

(b)  f/g  <  0  (for  softening)  ,  (43) 

(c)  f/g  =  0  (for  perfectly  plastic) 

A 

We  emphasize  that  a  condition  of  loading,  i.e.,  g  =  0  and  g>0,  is  always  pre¬ 
supposed  in  the  definitions  (43).  It  is  worth  observing  from  (24b)^  that  once 

PKL’  ^KL*  g  and  ^1  are  sPecified»  then  the  strain-hardening  response  is  also  known. 

We  now  provide  a  geometrical  interpretation  of  the  definitions  (43).  We 

A 

recall  that  during  loading,  since  g  =  0,  g>0  and  g  =  0,  the  strain  trajectory  Ce 

is  intersecting  the  yield  surface  3£  and  locally  pushing  it  outwards.  Since 

m  • 

g  =  0  and  g  =  0  it  follows  from  (17)  and  (19)  that  f  =  0  and  f  =  0  also,  and  the 

corresponding  stress  trajectory  C  is  intersecting  the  yield  surface  39  in 
stress  space.  If  the  material  is  hardening,  (43a)  holds  and  the  stress 
trajectory  Cg  is  directed  outwards  and  is  pushing  the  surface  39  locally 
outwards.  But,  (43b)  holds  if  the  material  is  softening  and  the  stress 
trajectory  is  directed  inwards  and  is  pulling  the  surface  39  locally  inwards. 

In  perfectly  plastic  behavior  when  (43c)  holds,  the  stress  trajectory  con¬ 
tinues  to  lie  on  the  yield  surface  39  which  is  stationary. 

Thus  while  during  loading  the  stress  trajectory  Cg  is  always  pushing  the 

*Since  g  is  always  positive  in  (43),  we  could  use  oi^Ly  f  in  providing  the 
above  definitions.  But  the  use  of  the  quotient  r/g,  which  is  rate- 
independent,  is  preferable  in  general.  For  certain  purposes,  however,  it 
is  useful  to  employ  only  f  as  in  (58)3  and  (59)  of  Section  4. 


yield  surface  BC  in  strain  space  locally  outwards,  the  corresponding  yield 
surface  B8  in  stress  space  may  be  moving  concurrently  outwards,  inwards  or  may  be 
stationary  depending  on  the  type  of  strain-hardening  response  being  exhibited.  The 
actual  occurrence  of  such  behavior  has  been  indicated  in  Section  1  with  reference 
to  the  simple  tension  test.  The  usual  stress  space  formulation  of  plasticity 
theory  introduces  a  priori  loading  criteria  in  stress  space  and  stipulates  that 
during  loading  the  yield  surface  in  stress  space  can  never  move  inwards. 

Viewed  in  the  context  of  the  present  development,  the  usual  stress  space  formula¬ 
tion  of  plasticity  is  seen  to  include  only  a  hardening  type  response  and  to 

* 

exclude  softening  and  perfectly  plastic  responses  .  Figure  2  illustrates  the 
three  types  of  material  behavior  discussed  above. 

The  definitions  for  hardening,  softening  and  perfectly  plastic  behavior  intro¬ 
duced  in  (43)  require  the  use  of  yield  surfaces  both  in  strain  space  and  stress 
space.  However,  it  may  be  noted  that  our  terminology  for  softening  and  hardening 
seems  to  be  consistent  with  the  geometrical  sense  of  these  terms  employed  in  a 
stress  space  formulation  by  Melman  and  Drucker  [8];  see  Fig.  5  of  their  paper. 
Also,  Prager  [9]  employs  the  terms  hard  and  soft  with  reference  to  material 
behavior,  but  his  sense  of  these  terms  differs  from  ours:  In  [9]>  a  hard 
material  is  one  whose  stress -strain  curve  always  lies  above  a  given  straight 
line  (representing  linear  elastic  response)  with  the  deviation  from  linear 
behavior  increasing  for  larger  deformation;  a  soft  material  is  one  whose  stress - 
strain  curve  always  lies  below  the  straight  line  with  the  deviation  increasing 
for  larger  deformation. 

In  what  follows,  we  frequently  need  to  refer  to  a  set  of  conditions  which 
must  be  satisfied  by  various  functions  and  material  coefficients,  and  which 
arise  from  characterization  of  strain  hardening  response.  To  avoid  undue 

*In  the  context  of  the  present  paper,  it  is. not  possible  to  formulate  loading 
criteria  in  stress  space  using  only  f  and  f. 


19. 


repetition  we  denote  this  set  of  conditions  by  H  and  write 


H 


{>  0  if  and  only  if  the  material  is  hardening  , 

<  0  if  and  only  if  the  material  is  softening  , 

=  0  if  and  only  if  the  material  is  exhibiting 

perfectly  plastic  behavior 


(a) 

(D> 

(c) 


(44) 


Returning  to  the  definitions  (43)  and  recalling  (24b )^  and  (12)25  it  is  seen 
that 


satisfies  conditions 


H  . 


(45) 


It  is  worth  mentioning  that  the  usual  treatment  of  an  elastic-perfectly  plastic 

material  (see,  for  example,  [2,  Sec.  9])  in  stress  space  requires  the  use  of  a 

yield  condition  of  the  form  f (s,_  )  =  const,  and  the  quantity  on  the  left-hand 

side  of  (45)  indeed  vanishes  identically  in  this  case. 

With  the  use  of  the  definitions  (43a, b),  we  now  obtain  an  expression  for 

the  rate  of  plastic  strain  which  is  valid  in  regions  of  hardening  and  softening 

behavior  only.  Thus,  by  (lOd),  (12),  (24b)1  and  (43a),  in  a  region  of  hardening 
•  n  A 

e£T  can  be  related  to  f  through  the  expression 


(46) 


A  •  p 

with  (43a)  and  (45a)  holding  ,  while  in  a  region  of  softening  e^  is  again 
given  by  (46)  but  now  with  (43b)  and  (45b)  holding;  in  both  cases,  the  sign  of 


t 


The  equation  number  (45a)  refers  to  (45)  along  with  part  (a)  of  condition  H. 


TABLE  1 

Summary  of  loading  criteria  in  strain  space  and 
associated  conditions  in  stress  space 


Elastic 

g<  0 

implies 

f  <  0 

A 

A 

Unloading 

g  =  0,  g  <  0 

implies 

f  =  0,  f  <  0 

Neutral 

Loading 

o 

ii 

<  hO 

•S 

o 

II 

implies 

A 

f  =  0,  f  =  0 

the  coefficient  of  p,_  in  (46)  is  positive.  For  perfectly  plastic  behavior,  it 

IUj 

is  clear  from  (lOd),  (24b).,  (43c)  and  (45c)  that  e£_  cannot  be  expressed  as  a 

1  All 

A  , 

product  involving  f  and  must  be  calculated  from  (10d).  For  convenience,  a 
summary  of  the  relationships  between  the  loading  criteria  in  strain  space  and 
the  associated  conditions  in  stress  space  is  provided  in  Table  1. 

In  the  remainder  of  this  section,  we  discuss  some  special  cases  of  the 
foregoing  results  which  are  of  particular  interest  in  view  of  their  simplicity. 
The  first  two  of  these  (see  cases  (a)  and  (b)  below)  examine  the  consequences 
on  strain  hardening  behavior  of  certain  restrictions  on  the  stress  response 
functions  in  (5)^  and  s^  in  (34).  The  third  (see  case  (c)  below)  pertains 
to  a  limiting  behavior  of  strain  hardening  response,  i.e.,  saturation  hardening 
and  softening. 

(a)  Consider  the  special  case  of  (5)  for  which  the  stress  response  is 


independent  of  its  last  two  arguments,  i.e., 


A 

as 


MN 


A 

as 


=  0 


MN 


ae 


3K 


=  0 


(47) 


KL 


Then,  by  (23)  we  have 


_26 _ K. 


ae 


KL 


ae' 


aa  _  af 

*  dK  3(C 


KL 


and  hence  by  (24a)  or  (22) 


A.A 

f/g  =  1 


(48) 


Recalling  the  definitions  (43),  it  is  clear  that  a  material  for  which  (47)..  _ 
hold  can  never  exhibit  softening  or  perfectly  plastic  behavior.  If  conditions 
(47)  are  satisfied  and  if  ag/ae^^O,  it  follows  from  (27)  and  (30)  that 

V*  =  0  ,  v  =  0  .  (49) 

Conversely,  if  (49)^  is  satisfied  so  is  (49)g  and  then  (48)  holds  by  virtue  of  (32). 
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(b)  Consider  the  special  case  of  (34)  for  which  the  stress  response  is 
independent  of  its  last  two  arguments.  In  this  case,  the  results  (36)  to  (42) 
hold.  It  then  follows  from  (37),  (40),  (43)  and  (12 )g  that 

T  satisfies  conditions  H  .  (50) 

With  the  use  of  (36)^  and  (40),  in  a  region  of  hardening  or  softening  (46) 
becomes 


ep  =  1  A  0 
“  r*KL 


(51a) 


while  it  follows  from  (50c),  (42),  (36^  and  (lOd)  that  in  a  region  of  perfectly 

plastic  behavior 


=  0  * 


'KL  A  as 


(51b) 


KL 


(c)  Caulk  and  Naghdi  [5]  have  previously  introduced  a  definition  of 
saturation  hardening  in  connection  with  their  discussion  of  hardening  response 
in  cyclic  loading  of  metallic  materials  (see  Eq.  (19)  in  [5])-  In  view  of  the 
definitions  (43),  it  is  of  interest  to  reexamine  here  the  notion  of  saturation 
hardening.  Thus,  for  our  present  purpose,  an  elastic-plastic  material  is  said 

to  exhibit  saturation  hardening  along  a  strain  trajectory  Cg  (or  a  stress 

.  & 

trajectory  Cg)  if  and  only  if  there  exists  a  constant  K^  such  that* 

lim  f/g  =  \  >  0  (g  =  0,  g > 0)  .  (52a) 

t-*® 


Similarly  a  material  exhibits  saturation  softening  along  a  strain  trajectory 

C  if  and  only  if  there  exists  a  constant  K  such  that 
e  s 

lim  f/g  =  K  <  0  (g  =  0,  g>0)  .  (52b) 

t-.®  s 

A 


In  the  definitions  (52a, b)  we  have  excluded  for  convenience  the  equality  sign. 
If  the  limit  of  the  left-hand  sides  of  (52a, b)  is  zero,  we  say  that  the 
material  saturates  to  a  perfectly  plastic  behavior. 


4.  Strain  Hardening  Response  for  Special  Constitutive  Equations. 

We  consider  now  in  some  detail  the  nature  of  the  hardening  response  in 

small  deformation  of  metals  whose  behavior  is  characterized  by  a  simple  set 

of  constitutive  equations  appropriate  for  elastic-plastic  materials  which  are 

homogeneous  and  initially  isotropic  in  their  reference  state.  First,  we 

6  P 

recall  that  the  infinitesimal  elastic  strain  tensor  is  defined  by  eKjj  =  eKL~eK 

P  G 

and  note  that  with  efL  = 0  in  the  reference  configuration,  e  =0  also  there. 

KL  KJj 

It  is  convenient  to  utilize  a  standard  decomposition  for  second  order 
tensors.  Thus,  for  example,  in  the  case  of  the  stress  tensor,  we  have 


SKL  S6KL+TKL  ’ 


where  s6,_  is  the  spherical  part  of  s„T,  t 

KJ_i  Mj  rJj 

part  of  s  and  s  is  the  mean  normal  stress. 

KJj 

eKL’eKL’eHL  int°  sPherical  Parts  ^KL’^KL’® 


■  1  c 

3  SKK  * 

is  the  deviatoric  (traceless) 

In  a  similar  manner  we  decompose 
^6^  and  deviatoric  parts  y^. 


YKL,YKL' 

Let  the  stress  response  function  in  (5)^  be  specified  by  generalized 
Hooke's  law,  namely 


tkl  *  '  s  ■  3k?!  >  <53> 

and  the  coefficient  function  for  the  rate  of  work-hardening  response  in 
(9)  in  the  form  [10] 


CKL  ■  eTKL  +  WtKL 


(54) 


where  ^  is  the  shear  modulus,  k  the  bulk  modulus  and  j3  and  0  are  constants. 

With  the  use  of  the  decompositions  just  noted,  the  loading  functions  f(\j)  and 
g(U)  can  be  written  as  different  functions  f(TMN>s and  g(yMN>e,y^N,eP,ic). 
In  this  section,  we  restrict  attention  to  special  loading  functions  of  the  form 
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W  Jl  «H» 


■*« 


f(u)  =  =  Tj^Tj^ +  3^-10  , 

S(U)  =  6(YMI,.e»'vSI,»eP.K)  =  ^2(vKL-Y^L)(vKL-y^L)  +27^2(e-?>)2  -  K  , 

where  ji  is  a  constant  and  where  (17)  and  (53)  have  been  used^.  Utilizing 
formulas  of  the  type"*' 


(55) 


af 


as 


MN 


af  i/  af  af n . 

3tmn'3  SocV*1* 


and  recalling  (20)  and  (ll),  it  can  be  easily  shown  that 

isTT  =  2^TMN+ 'lrS6MN^  *  f  =  2^tmntmn+3,|,s  ’ 

MN 

8  =  2(2^WMN  +  9*k?^  =  ^  +  2  (2M-tmmYmw  +  9^‘®f>) 


MN  MN 


(56) 


(57) 


and  the  expressions  for  ag/de^  ar.d  df/de^  may  be  obtained  similarly.  We 

A  •  • 

recall  that  during  loading  g  is  positive  while  g=f  =  g  =  f  =  0.  Keeping  this 
in  mind,  it  follows  from  (55 )1  and  (57)2  that  during  loading 

tkltkl+3,|,s  "  *  =  0  >  2(tki/kl  +  3* s  "  *  =  0  »  f  =  k  (58) 

and  hence  by  the  definitions  (43), 

K  and  (tmhtmn + 31* s s )  both  satisfy  conditions  H  .  (59) 

Clearly  for  the  special  constitutive  equations  used  in  this  section,  in  view  of 
(58)^  and  (59) >  the  strain  hardening  behavior  may  be  characterized  by  K.  Further¬ 
more,  frcm  (58)^  and  (10b, c)  during  neutral  loading  it  is  necessary  that  k=0  and 
during  unloading  it  is  necessary  that  K<0.  In  this  connection,  recall  Table  1 
and  the  discussion  following  (21). 


^The  loading  function  (55)^  does  not  depend  explicitly  on  plastic  strain,  but 
includes  a  dependency  on  mean  normal  stress.  When  i|r  =  0  and  K  =  const.,  (55)^ 
reduces  to  the  usual  von  Mises  yield  function.  A  loading  function  of  the  type 
(55)^  was  previously  employed  by  Green  and  Naghdi  [10]. 

^It  is  understood  that  in  line  with  the  summation  convention,  our  notation 
a^TKK  in  (56)  stands  for  the  sum  af/dT^  +  + 
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The  stress  response  (53)  nay  be  regarded  as  a  special  case  of  that  in 
(3*0  with  the  last  two  arguments  absent;  and,  in  addition,  the  symmetry  con¬ 
ditions  mentioned  following  (34)  are  satisfied  by  (53)*  Hence,  in  addition  to 
(36)  to  (42)  the  special  results  obtained  at  the  end  of  Section  3  [see  case  (b) 
following  Eq.  (49)]  remain  valid  here.  Thus,  using  (54)  and  (55)>  from  (33)> 

(4l)  and  formulas  of  the  type  (56)  and  (57)^»  we  obtain 

A  =  4(2htkltkl  +  9t2kSe)  ,  r  =  2(pTKLTKL  +  3^i2)  .  (60) 

With  the  use  of  (60),  \y  and  f/g  may  be  obtained  at  once  from  (42).  Also, 
remembering  (50),  we  observe  that  in  this  case  the  right-hand  side  of  (60)^ 
provides  a  rate-independent  characterization  of  strain  hardening.  Constitu¬ 
tive  equations  for  the  rate  of  plastic  strain  or  equivalently  for  and 
simplify  and  may  now  be  obtained  from  (51a)  in  a  region  of  hardening 
or  softening  and  from  (51b)  in  a  region  of  perfectly  plastic  behavior. 

Since  our  development  in  Sections  2  and  3  began  with  the  strain  space 
(rather  than  the  stress  space)  formulation  as  primary  and  since  the  quotient  f/g 
is  used  to  define  strain  hardening,  it  is  desirable  to  examine  the  predictions 
of  various  theoretical  results  in  the  case  of  the  familiar  one -dimensional 
tension  test.  To  this  end,  consider  a  homogeneous  deformation  sustained 
by  a  uniaxial  tension  s^  =  s  =s(^)  along  the  X^-axis.  Then,  using  a  matrix 
representation  for  t^,  we  have 

2  0  0 
0-10 
0  0-1 

where  for  brevity  we  have  introduced  the  constant  matrix  ))b^J|.  Assuming 

that  initial  yield  occurs  at  a  value  s  of  s  and  a  value  K  >0  of  K>  the 

00 

solution  can  be  obtained  in  a  straightforward  manner.  We  omit  details,  but 
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record  here  some  of  the  results  of  interest  : 


(C  =  js2(2+*)  ,  Ko  =  |-s2(2+t)  ,  2  +  *>0  ,  6q>0  ,  (6 

s  >  0  ,  K  >  0  ,  e^  >  0  when  g  =  0,  g  >  0  ,  ( 6; 

Both  s  and  (2p+\|f0)  satisfy  condition  H  .  (6 


We  postpone  a  discussion  of  perfectly  plastic  behavior  until  later  in  this 
section  but  consider  further  calculations  for  the  other  two  types  of  behavior: 
In  a  region  of  hardening  or  softening,  the  elastic  and  plastic  strains  are 


-  iwy  - 1  m  > 

?  =  ^  .  11411 -7^  IW  . 

9k  6p, 


de 


11  dYll 
ds"~  and  ”ds“ 


>  0  if  and  only  if  the  material  is  hardening 


<  0  if  and  only  if  the  material  is  softening 


*  *  *  * 

where  the  constants  v  ,E  ,pi  and  k  are  defined  by 


* 


v 


* 

k  = 


_  hi 

2+t 


y 


* 

E  = 


2(ln>*) 


3(l-2v  ) 


ill 

9+ 


3(2B+*d) 

2(2+t)2 

*?«*  • 

E*  U/0) 


(6 


(6 


and 


"it  is  clear  from  (58)3,  (10),  (62a  )i  3  and  (62b)x  that  during  neutral  loading, 
it  is  necessary  that  3  s  =  0  and  duri&g  unloading  it  is  necessary  that  s<0. 
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,'****  V*  V-.'*  -v 


*  # 

E  and  ^  satisfy  conditions  H 


(63c) 


The  constants  in  (63d)  have  been  defined  analogously  to  the  corresponding 


E  E 

constants  in  linear  elasticity,  e.g.,  M.  =  ^'(l+v)'  »  k  =  3(i-2v)  *  where  v  is 


Poisson's  ratio.  In  the  special  case  that  i|r  =  0,  v  and  the  expressions 


*  * 


for  E  simplify  while  k  as  ^-*0.  Continuing  our  discussion  of  haidening 


and  softening  behavior,  it  can  be  shown  that  when  \|r  ^  0  (see  Appendix  A  for 

A  ,A 


details)  the  quotient  f/g  can  be  written  as 

_ 2  +  t 


A 

f 

A 

g 


3  ~"*'i  ds  11  11  ds  11  ‘  *  ds  '  ds  '  nuKL’ 

* 

where  tr  stands  for  the  trace  operator.  In  the  special  case  when  v  =v, 
(64)  reduces  to  (see  Appendix  A  for  details) 


I-  ce g-]-1  =  [l+^r1  , 


and  by  (43a, b) 


de 


,  de 

and  (1 +  -=-£) 
ds  de 

e 


>  0  if  and  only  if  the  material  is  hardening  , 
<  0  if  and  only  if  the  material  is  softening  , 


where  as  in  Section  1,  we  have  again  used  the  notation  e = e, , ,  e  = ef , , 

11  e  11’ 


e  =  e?\  . 
P  11 


it  is  clear  that  during  loading  e^  is  strictly  increasing  with  time  by 


virtue  of  (62b)^.  Moreover,  according  to  (62c)  the  time  rate  of  stress 


may  be  used  to  characterize  strain  hardening  behavior  in  uniaxial  tension 
and  a  characterization  of  the  same  behavior  is  provided  by  the  combination 


(64) 


(65a) 


(65b) 


Before  closing  this  section,  it  is  desirable  to  elaborate  briefly  on 
some  features  of  the  foregoing  results  for  uniaxial  tension,  which  have  been 
obtained  with  the  use  of  a  special  set  of  constitutive  equations.  With 
reference  to  all  three  types  of  strain  hardening  response  defined  in  (43), 
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(20+^0)  of  the  constitutive  coefficients.  While  the  elastic  moduli  E,n  are 
always  positive,  it  follows  from  (63e)  that  the  constants  E  ,p,  are  positive 

in  a  region  of  hardening  and  negative  in  a  region  of  softening.  In  the  special 

*  /r  \  A ,A 

case  of  v  =v,  it  is  clear  frcm  (65a)  that  the  quotient  f/g  can  be  expressed  in 

terms  of  quantities  (2)  to  (4)  and  indeed  (65b)  corresponds  to  the  behavior 

summarized  in  (4)  for  uniaxial  tension^.  Furthermore,  with  f  =  0  in  (55)-^  the 

a 

plastic  volume  change  or  equivalently  vanishes  also.  The  strain-hardening 

response  is  then  characterized  by  0,  in  view  of  (62c).  Also,  in  a  region  of 

A,A  .  . 

hardening  or  softening  the  quotient  f/g  reduces  to  (see  Appendix  A  for  details) 


g 


e  e 

where  we  have  put  Y=Y11  and  Y  =  Y-y_» 

The  significance  of  the  strain  space  formulation  in  the  case  of  elastic- 

r  1  A, A 

perfectly  plastic  materials  was  pointed  out  in  [1J.  Since  the  quotient  f/g 

is  used  here  to  define  various  types  of  hardening  response,  it  is  desirable  to 

indicate  the  reduction  of  the  present  development  to  the  usual  perfectly  plastic 

.  A  . 

behavior  in  uniaxial  tension.  First,  we  observe  that  during  loading  (g  =  0,  g>0) 
for  perfectly  plastic  behavior  f/g=T=0  by  (43c)  and  (42).  It  then  follows  that 
K  =  0,  tc=K  ,  s  =  s  by  (58)  and  (62a).  0  and  that  el?..,  although  indeterminate*, 
is  strictly  increasing  with  time  in  view  of  (626)^.  Thus,  in  the  context  of  the 
present  paper,  the  uniaxial  stress -strain  curve  for  elastic-perfectly  plastic 
behavior  consists  of  a  linearly  elastic  portion  followed  by  a  horizontal  portion 

and  as  time  progresses  the  locus  of  e^  moves  outward  along  the  abscissa  of  the 

^Recall  that  the  special  constitutive  equations  employed  in  this  section  are 
not  sufficiently  general  to  predict  all  details  of  the  stress -strain  curve 
in  Fig.  1.  Indeed,  different  choices  of  the  combination  (20+$0)  of  the 
coefficients  (appropriate  for  different  materials)  yield  stress-strain 
curves  consisting  of  straight  line  segments  whose  slopes  correspond  to  the 
rising  or  falling  portions  of  the  curve  in  Fig.  1. 

*The  indeterminacy  of  e-Q  stems  from  the  fact  that  (51b)  in  this  case  reduces 
to  an  identity. 


s-e  curve.  This  is  in  agreement  with  the  usual  characterization  of 

perfectly  plastic  behavior  in  uniaxial  tension.  We  also  note  that 

an  examination  of  the  solution  given  by  (63a, b)  and  (63e) 

easily  reveals  that  hardening  (softening)  is  represented  in  a  stress-strain 

diagram  by  a  straight  line  which  lies  above  (below)  the  horizontal  perfectly 

'if 

plastic  line.  Indeed,  since  e^  =  s/E  +  (s-sq)/E  ,  then  de^/ds  =  l/E  +  l/E 
and  by  standard  results  for  inequalities  it  follows  from  (63e)  that 

^ell  1  1 

oo  >  ■  >  max{— ,  — }  if  the  material  is  hardening  , 

u"  E 

1  ^ell  1 

—  >  — t —  >  — i  >  -oo  if  the  material  is  softening 
E  ds  E* 

Moreover , 


dell  1 

oo  >  — —  >  —  implies  that  the  material  is  hardening  , 
ds  ll 


dell  1 

-oo  <  •  <  -  implies  that  the  material  is  softening  . 


5.  Saturation  hardening 

As  in  Section  4  we  again  restrict  attention  to  small  deformations  of 
elastic-plastic  materials ,  which  are  homogeneous  and  initially  isotropic  in 
their  reference  configuration.  We  also  assume  that  there  is  no  plastic 
volume  change  so  that  ?  =  0  in  the  notation  of  Section  4.  For  a  fairly  large 
class  of  metallic  materials,  it  is  well  known  that  the  stress -strain  curves 
of  uniaxial  cyclic  loading  attain  --  after  several  cycles  — 
saturation  hardening.  The  purpose  of  this  section  is  to  indicate  how  the 
development  of  Sections  2  and  3  can  be  used  to  characterize  a  hardening 
response  that  includes  saturation  behavior  and  to  compare  the  results  with 
those  of  Caulk  and  Naghdi  [53* 

Starting  with  a  fairly  general  discussion  of  loading  functions  contained 
in  the  paper  of  Green  and  Naghdi  [2],  for  initially  isotropic  materials  Caulk 
and  Naghdi  [5]  derived  a  loading  function  in  the  form  (see  [5,  Eqs.  (40^  and 
(56)l3) 

f(\j)  =  f( tMN’^MN’*)  =  tkltkl  ”  “tklykl  +  ctvklykl  "  K  * 

(67) 

g(u)  =  =  4u  (Y^-Y^)(y^-Y^)  -2o^(v^-V^L)Y^L  +  cry^LY^L- k  > 

where  a  and  a  are  constants  and  where  (53)  has  been  used  in  writing  (67 )g.  It 

should  be  noted  that  the  loading  functions  (67).  „  depend  explicitly  on  but 

not  on  the  mean  normal  stress  s.  Here  we  also  adopt  (67),  0  but,  instead  of  the 

1,2 

hardening  response  assumed  in  [53  >  we  specify  the  coefficient  function  Cj^  in  (9)  by 

CKL  =  9(<)tKL  +  ^(k)y?L  ’  (68) 

which  is  different  from  that  used  in  Section  4.  The  constitutive 

A  A 

assumption  for  in  C 5 3  is  similar  to  (68)  but  with  p(<)  and  T](k)  specified 
by 
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(69) 


POO  = 


K-tt  a  tC-K 

S  A  .  s 

U  P  >  Tj( K J  =  u  T|  > 

OS  O  b 


where  0  and  T|  are  constants,  tco  is  the  value  of  k  at  initial  yield  and 

K  is  the  saturation  value  of  K.  Since  the  stress  response  (53)  is  used  in 
s 

this  section,  in  addition  to  (36)  to  (42),  all  the  results  stated  under 
case  (b)  at  the  end  of  Section  3  are  also  valid  here. 

To  facilitate  the  discussion  that  follows  and  for  later  reference,  we 
record  the  expressions 

aT“  =  2tMN  "  “^MN  ^  0  »  =  2n  /  °  ,  (g  =  0)  , 

dMN  ^  aeMN  ^MN 

f  =  ^2TMN-arYMN^TMN  * 

g  =  =  f  +  2n(2TMW-Qfy^M)v^w  , 


MN  TMN'  'MN 


which  have  been  obtained  with  the  use  of  formulas  of  the  type  (56)  along  with 
(20),  (11),  (36)1  and  (37).  With  the  help  of  (70)  and  recalling  the  definitions 
(33)  and  (4l),  A  and  r  are  given  by 

A  =  2,i  as  as  =  2^^2tkl-q^^2tkl'“ykl^  >  0  * 

K-L  Mj  ^  rj2 

r  =  (2TKL-Q'YKL)^(0f+^(,t))TKL"  (2«y-T|(K))v|L}  , 


Thus,  based  on  the  constitutive  equations  assumed  in  this  section, 

A  is  2 (j,  times  the  square  of  the  magnitude  of  the  normal  to  the  yield 

—  * 

surface  a§  in  stress  space.  Having  obtained  the  results  (71 W  can  be 
calculated  from  (42  )^  and  it  then  follows  from  (32)  that  the  quotient  f/g  must 


satisfy  the  inequality 


A.A 

f/g  <  1  , 


which  limits  the  extent  of  the  hardening  behavior.  The  restriction  (72),  in 


1  ,rj.  V*  ‘4 

^  /. it n  j-, - .A'i latur  i  r*^ 


•W^linCr  r  * 


turn,  places  an  upper  bound  of  unity  on  the  value  of  the  saturation  constant 
in  (52a)  so  that 

0<KhSl  .  (73) 

Expressions  for  can  now  be  easily  calculated  from  (51a)  in  a  region  of 
hardening  or  softening  and  from  (51b)  in  a  region  of  perfectly  plastic 
behavior. 

Given  the  constitutive  assumptions  employed  in  this  section,  the  results 
(71),  0  and  the  restrictions  (72),  (73)  are  valid  for  any  small  elastic- 
plastic  deformations.  In  the  rest  of  this  section,  however,  we  again  confine 
attention  to  a  homogeneous  deformation  sustained  by  uniaxial  tension  (6l). 

Since  plastic  volume  change  e^=0,  e6  =  e  is  given  by  (63a)^.  Again,  as  in 
(l)^,  for  convenience  we  use  the  notation  e  =  e^,  ep  =  e,)  1 ,  e^e.^  and 
write 

IWklW  ~  *  ep^bKL^  * 

where  the  constant  matrix  Hb^JI  is  defined  by  (6l)^.  Also,  from  (70)^^  4> 

(61)  and  (1)2*  we  deduce  that 

f  =  (js-orep)s  ,  g=  (j  s  -  ofep){Ee  +  (3^-E)ep)  ,  j  s  -  orep  /  0  (g  =  0)  .  (74) 

At  initial  yield  ep  =  0  and  kq  =  |  s^>0  by  virtue  of  (74)^  and  (67 )1. 

Hence,  on  the  yield  surface  (g  =  0),  |  s  -  crep  must  be  positive.  From  this 
last  result,  along  with  (30)  and  (36 )g»  we  have  ep>0  during  loading  and 
therefore  ep  is  strictly  increasing  with  time.  Further,  from  the  definition 
(43),  and  the  positivity  of  the  coefficient  of  s  in  (74)x,  it  follows  that  s 

* 

must  satisfy  the  conditions  in  (44).  The  above  results  may  be  summarized  as  follows 

*The  inequality  (75 )i,  together  with  (74K  and  (10),  imply  the  following: 

During  neutral  loading  it  is  necessary  that  s  =  0,  while  during  unloading  it 
is  necessary  that  s<0. 


(75) 


j  s  -  crep  >  0  ,  ep  >  0  , 


s  satisfies  conditions  H 


While  (75)2  holds  during  all  three  types  of  strain-hardening  behavior  defined 
in  (43),  it  follows  at  once  from  (75 and  (l)2  that  eg  also  satisfies  con¬ 
ditions  H. 

For  uniaxial  tension  under  discussion,  the  quantities  A  and  r  in  (71^  2 
reduce  to 

A  =  3n(j  s  -  aep)2  >0  ,  r  =  |  ( j  b  -  arep)F  , 

(1 

F  =  |  (flf+p(K))8  -  (2<j-fi(ic))ep  , 

where  for  later  convenience  we  have  introduced  the  quantity  T  defined  by  (76)^* 
Further,  from  (50),  ( 76 and  (75 )1  follows  the  result 


T  satisfies  conditions  H  . 


Also,  the  expression  for  plastic  strain  rate  ep  in  a  region  of  hardening  or 

softening  can  be  written  as 

|(|s-ae).  |E(|s-ofe)  . 

e  =  J  A  _ - 2-  s  =  J  ■  ■?  ,  ■£ - e  ,  (7 

p  r  r  +  |  E(|  s  -orep) 

which  is  similar  in  form  to  that  obtained  in  [5]  and  where  the  relation 

s  =  E(e-ep)  has  been  used  in  deriving  (78)2*  In  fact,  if  the  coefficient 

A  A  —  \ 

functions  0  and  T|  which  occur  in  r  are  specialized  to  those  given  by  (69)  >  then 

(78)  reduces  to  that  in  [5,  Eq.  (80)]. 

The  result  (75)2  enables  us  to  calculate  the  slopes  de/ds,  dep/ds 

explicitly  as  functions  of  s,ep,tC.  Thus,  with  the  use  of  (l),  (78^  and  chain 

rule  of  differentiation,  in  a  region  of  hardening  or  softening  we  have 


•  » 


r 


si  =  i+!!e  aep  H'-V  (, 

ds  E  ds  *  ds  jr  ' 

It  follows  from  (79),  (75)L  and  (77)  that 
de  de  1 

®  >  >0  ,  ®  >  ~  >  max{^,  >  0  if  the  material  is  hardening  , 

de  ^  .  de  * 

-®  <  -j-£  <o  ,  —  >  —  >  if  the  material  is  softening  . 

Moreover , 

de  1  116 

®  >  >  |j  (or  equivalently  ®  >  >0)  implies  hardening  , 

de  1  de  ^ 

-®  <  —  <  —  (or  equivalently  -®  <  -£■  <  0)  implies  softening  . 

Since  de  /ds  =  ^  which  is  always  positive,  we  may  write  de  /de  =  E  de  /ds, 
de/dee  =  E  de/ds  and  then  obtain  explicit  expressions  for  these  derivatives 
from  (79)*  It  is  evident  that  conditions  of  the  type  indicated  in  (80)  for 
dep/ds  also  hold  for  de^/de^  It  follows  from  (42),  (76)  and  (79 )g  that  in  a 
region  of  hardening  or  softening 

f  -i  *^ii  ^Tii  -i  -i 

1>£-  J  •  <f 

g 

In  view  of  (80a)  and  (43),  (8l)  implies  that  de^/ds  <-l/3n  in  a  region  of 
softening.  It  is  clear  from  (8l)^  and  (l)  that  a  knowledge  of  ^,E  and  the 
slope  de/ds  suffices  to  determine  f/g.  If  the  material  saturates  to  perfectly 
plastic  behavior,  the  left-hand  side  of  (81)^,  i.e.,  £/g  must  tend  to  zero  and 
hence  in  this  case  de^/ds  must  become  unbounded. 

We  now  turn  to  a  brief  discussion  of  saturation  hardening  usually 
observed  under  uniaxial  cyclic  loading.  Recalling  the  definitions  (52a, b), 
from  (8l)  we  deduce  that  saturation  hardening  occurs  if  there  exists  a  constant 
K.  such  that 


(80a) 


(80b) 


■h 

■u 


=  ,  (OO^Sl)  . 


t  -*® 


"t-»® 


(82) 


In  order  to  exploit  the  implications  of  (82),  we  first  observe  that  f  defined 


by  (76)3  can  be  rewritten  as 


r  =  |  (j  s  -  erep)(af+p(ic))  +(  2a  +  T\(k)  +|  p(tc))e{ 


(83) 


and  then  express  (79 )  in  the  form 


de  o  a  (a2 -2o  +  ^(ic) +|  P(ic))e  , 

-£=  [J  (a+J(ic))+ - “ — p - 2 - £]  1 


ds 


T<?s  •  ”p> 


(8M 


Consider  now  a  special  material  response  which  corresponds  to  the  vanishing  of 
the  numerator  of  the  second  term  in  the  square  brackets  in  (84),  i.e.. 


%--2o  +  l}(lc)  +|  pU)  =  0  , 


(85) 


which  has  the  same  form  as  a  particular  case  discussed  in  [5].  From  (75)1> 


(77)-,,  (83)  to  (85)  and  (8l),  it  can  be  readily  concluded  that 


a+p(k)  satisfies  conditions  H  , 

4/3 


^!e  . 

u«> 


>  f-  [i+- ^-r1  <1  . 

<*+P(tc) 


(86) 


A 

g 


Also,  in  view  of  (*'38 )  and  (86)g,  in  a  region  of  softening: 


0  >  a  +  p(ic)  >-  4p,  . 


If  saturation  hardening  occurs  with  0<Kh<l,  then  from  (82 )^,  (86)  and  the 
condition  (85)  we  have 


A  VS,  A 

lim  P(k)  *  a  ’  11a  T\M  =--irri  +  2a  ,  (0<Kh<l)  , 

-  n  t-»«  n 


(87) 


t,-»® 


36. 


while  g(tc)  becomes  unbounded  for  K^*l. 

We  further  examine  saturation  hardening  by  adopting  the  special  coefficients 
(69)  subject  to  the  condition  (85).  When  saturation  is  assumed  to  occur,  the 
limit  of  the  coefficients  (69)  as  t-*®  is  zero  and  from  (85),  (82)  and  (86)^ 
we  obtain 

cr2  =  4o  ,  ag  +  2T|  =  0  ,  lim|f=|+^  ,  0<  =  [1  +  <1  ,  cr>0  ,(88 

—*00 

the  first  three  of  which  are  the  same  as  those  derived  in  [5,  Eqs.  (70)  and 

(86)]. 

By  way  of  illustration,  consider  the  304  stainless  steel  whose  behavior 

in  cyclic  tens ion-compress ion  is  discussed  in  [5,  Sec.  7].  As  in  [5]>  for  the 

304  stainless  steel,  we  take  the  values  of  E=123  GRa  and  de/ds  =  (3.85  GFfe)  1 

at  initial  yield  and  also  assume  the  value  v  =  0.3  for  Poisson's  ratio*.  With 

these  values,  the  expressions  (79)^  and  (8l)  predict  that  the  quotient  f/g  at 

initial  yield  is  approximately  equal  to  0.027 ♦  Again  using  the  above  values, 

as  well  as  or  =  1.5  (for  tension),  (88)..  gives  an  approximate  value  of  0.008  for 
A, A 

K^.  Thus,  f/g  decreases  from  a  value  of  0.027  at  initial  yield  to  a  value  of 
0.008  at  saturation.  It  is  clear  from  (82)  that  the  definition  of  saturation 
hardening  given  by  (52a)  implies  that  the  slopes  de/ds  or  ue^/ds  tend  to 
constant  limits  at  saturation .  In  this  connection,  it  should  be  noted  that 
when  g  and  T]  are  of  the  form  (69),  the  definition  of  saturation  hardening  used 
in  [5]  also  gives  constant  limiting  slopes. 

We  return  once  more  to  the  perfectly  plastic  case,  and  first  observe  that 
the  expression  for*  can  be  obtained  from  (51b)  with  the  use  of  (70)^  ^  and 
(71)-^.  In  view  of  (75)3*  s = 0  for  perfectly  plastic  behavior  and  s  retains  its 

*A  value  for  Poisson's  ratio  was  not  needed  in  the  calculation  given  in  [5» 

Sec.  7).  With  v  =  0.3  and  E=123  GRa,  p,  is  calculated  to  be  47.31  GRa. 

*In  fact,  in  the  case  of  uniaxial  tension,  the  resulting  expression  is  an 
identity. 


initial  yield  value  sq  and,  in  accordance  with  (75)g>  ep  *-s  strictly  increasing 
with  time  during  loading.  The  work-hardening  parameter  K  may  then  be  obtained 
as  a  function  of  ep  from  g  =  f  =  0  with  f  given  by  (67)^: 

K  =  |  s^  -  os  e  +  |  ael  .  (89) 

3o  o  p  2  p 

By  (76)  and  {Tic))  for  perfectly  plastic  behavior  it  is  necessary  that 

j  (a+3(tc))so  -  (2a-T)(tc))ep  =  0  (90) 

for  all  ep.  We  observe,  however,  that  in  view  of  (76)  and  (77c)  the  constant 
values 

3(k)  =  —  &  ,  tJ(k  )  =  2a  (91) 

are  sufficient  for  perfectly  plastic  behavior.  It  should  be  noted  that  the 
values  (91)  satisfy  the  condition  (85). 
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We  provide  here  details  of  the  calculations  leading  to  (64),  (65a)  and  (66), 

A. A 

and  also  record  alternative  useful  forms  of  the  quotient  f/g  associated  with  the 
constitutive  equations  of  Section  4.  From  (63),  in  a  region  of  hardening  or 
softening  we  have 


.  <*>“» 


dV,.T 

ihsrll 


dY 


KLii-l 

ds 


K.,,  4Y 


,e 

- IMI  *  ihs-Hl  Ihjrir1- 


(Al) 


(A2 ) 


de  de„j  i  de„.j 

11-dTll  ll-^rir1  =  +  ll-irli 


de 


KLi  -1 
ds  1 


I6klI!  +  7 


* . 

0  v  /v 


0 

0 


*  / 

0  V  /v 


,  (A3) 


In  a  region  of  hardening  or  softening  T  /  0  by  (50)  and  using  (42 )g  we  may 
write  f/g =  (1  -  jr)"1.  Then,  by  (60)  and  (61)  we  have: 

L  _  n+alVjai&l]-1 

£  ■  1  2p+*0  J 


=  [1  +  _E_  f +  l-i.-gy-If}]-1 

U  *  1  1-*V  l-2v 


3E 


A, A 


and  we  may  recall  that  1  +  v  >  0,  1  -  2v  >  0.  If  $/0,  then  f/g  may  also  be 
written  as 

£.  r,+  i-l  .  _ S+i -  . 

g  4/+3*V  2(1  +  -^)  +  4(1  +  4) 


(A4) 


(A5) 


The  result  (64)  follows  at  once  from  (Al),  (A2)  and  (A5).  Similarly,  if 
v* * v  (or  f  =  y“),  then  fl>om  (A3)  and  (A4)  we  obtain 


I-  ci-lr1-  tf  Mil. JMI^I  ll^irV1  • 

g  E 


(a6) 


The  result  (65a)  follows  from  (a6),  (63a)  and  (64).  With  the  use  of  (A3)  it  is 


de 


de' 


-  WT  1 

also  possible  to  write  (A6)  in  terms  of  H- r~ll  II  II  tut  we  do  not  record 


ds 


this  here. 

In  the  special  case  that  $  =  0,  we  note  that  by  (63d)2  ^  ai*d  (A2), 


1  +  4  J  =  itrt||dYKLM  '■dYKL"-1 


1  _  11  / dV  \  -1 

J  "  Hs  '  Hs  J 


ds  11  11  ds  11  J  ds  '  ds 
The  relations  (a4)x,  (A7),  (63a)g  and  (63b)2  lead  to  the  expression  (66). 


(A7) 


41. 


Fig.  1.  Idealized  stress -strain  diagram  for  a  typical  ductile  metal.  As  the 
points  1,2,3»4,5  of  the  stress-strain  curve  are  successively  traversed, 
the  locus  of  the  yield  point  on  the  e-axis  moves  outwards  through 
B1’B2,B3,B4  and  B5*  resPectively>  while  the  corresponding  locus  of  the 
yield  point  on  the  s-axis  first  moves  upwards  through  A^A^  to  A^>  and 
it  then  moves  downwards  through  and  A^.  All  unloading  curves  are 
drawn  parallel  to  the  linear  elastic  segment  1-0  and  hysteresis  is 
ignored. 

Fig.  2.  A  sketch  indicating  the  motion  of  yield  surfaces  in  stress  space  and 
strain  space.  During  loading  the  yield  surface  dC  in  strain  space  moves 
outwards  with  the  strain  trajectory  C  through  positions  such  as 
B^Bj^B^..  The  corresponding  yield  surface  dS  in  stress  space  moves  out¬ 
wards  through  positions  such  as  A^  and  A^  during  hardening  behavior,  is 
stationary  in  positions  of  the  type  A^  during  perfectly  plastic  behavior 
and  moves  inwards  through  positions  such  as  A^  and  A<_  during  softening 
behavior. 


